DETERMINAREA SOLUŢIILOR ECUAȚIILOR DIOFANTICE 
2069-2080 


CIRA OCTAVIAN 


REZUMAT. În acest articol ne propunem să determinăm soluţii 
pentru ecuaţiile diofantice 2069-2080. 


1. INTRODUCERE 


Funcţia care asociază fiecărui număr natural n pe cel mai mic nu- 
măr natural m care are proprietatea că m! este multiplu lui n a fost 
considerată de prima dată de Lucas, [3], în anul 1883. Alţi autori care 
au considerat această funcţie în lucrările lor sunt: Nueberg (1887) [4], 
Kempner (1918) [2]. Această funcţie a fost redescoperită de Smaran- 
dache, [5], în anul 1980. Funcţia este notată de Smarandache cu $ 
sau cu 7), iar în site-ul Wolfram Math World, [10], se notează cu pi. În 
acest articol am adoptat notația 7 a lui Smarandache din lucrarea [7]. 

În monografia [7], Smarandache prezintă o listă de ecuaţii diofan- 
tice referitoare la funcţia 7. Fie m,n € N* fixaţi şi x,y e N necunos- 
cutele. Ecuațiile diofantice 2069-2080 cu funcţia 7 sunt: 


(2069) n(m-z+n) =, 


(2070) n(m:z+n)=m+n-z, 
(2071) n(m-z+n)=az!, 
(2072) n(z”) =", 
(2073) n()” = n(z”), 
(2074) n(m- 2 + n) = (az), 
( 


(2075) n(z) + y = 2 + n(y), unde z şi y nu sunt prime. 

(2076) n(z) + n(y) = n(z + y), unde z şi y nu sunt prime gemene. În 
general funcţia 7 nu este o funcţie aditivă. 

(2077) n(z + y) = n(z) - n(y) . În general funcţia n nu este o funcţie 
exponențială. 

(2078) n(z - y) = n(2) : n(y) . În general funcţia 7 nu este o funcţie 
multiplicativă. 

(2079) n(m-z+n) =, 

(2080) n(z) - y = z:n(y), unde z şi y nu sunt prime. 
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2. ALGORITMUL KUMPNER 


Funcţia n : N* — N*, p(n) = m, unde m este cel mai mic număr 
natural care are proprietatea că n îl divide pe m!, (sau m! este mul- 
tiplu de n) este cunoscută în literatura de specialitate ca funcția lui 
Smarandache. Valorile funcţiei, pentru n = 1,2,...,18,sunt: 1, 2,3, 4, 
5,3, 7,46, 5, 11,4, 13, 7,5,6, 17,6 obţinute cu un algoritm ce rezultă 
din definiţia funcţiei n, după cum urmează: 


Programul 2.1. 


pir) = | form = la 
return m îf mod(m!,n) = 0 


Acest algoritm nu poate fi folosit pentru n > 19 pentru că nume- 
rele 19!, 20!, ...sunt numere cu mai mult de 17 cifre zecimale şi în 
modul clasic de calcul (fără o aritmetică a preciziilor arbitrare [9]) se 
vor genera erori de trunchiere. 

Kempner, [2], a dat un algoritm pentru a calcula pe p(n) folosind 
factorizarea clasică py" - p3? -: :p$* cu numere prime a lui n, şi baza 
de numerație generalizată (a;)y,, pentru i = 1,s. Soluţii parţiale 
pentru algoritmul de calcul a lui p(n) s-au dat mai devreme de Lucas 
şi Neuberg, [8] . 

Prezentăm algoritmul lui Kempner, [2], de calcul al funcţiei lui 
Smarandache 7. Pentru început definim şirul recursiv aj = p-a;+l 
cu j = 1,2,... şi a, = 1, unde p este un număr prim. Acest şir de 
numere constituie baza de numerație generalizată a lui p. Deoarece 
az =p+l,a3 =p?+p+1,...se poate demonstra prin inducţie că 
aj; =1l+p+...+pi1=(p'—1)/(p- 1) pentru Vj > 2. 

Valoarea lui v, astfel încât a, < a < au, este dat de formula 


(2.1) v = log, (1+a(p-—1))], 


unde |:| este funcţia parte întreagă inferioară. Cu ajutorul algoritmului 
lui Euclid putem determina şirurile unice K; şi r; după cum urmează 


(2.2) a = ka tru, 

(2.3) To > Ko Opal tul 

(2.4) 

(2.5) zii 3 Pearl 20 gat iat) leu fă 3 


(2.6) Tu=(A—1) > FA" OA: 
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adică, până când restul r,_x = 0. La fiecare pas i; este partea întreagă 
a raportului r;/a; şi r; este restul împărţirii. De exemplu pentru pri- 
mul pas avem k, = |a/a,] şir, = a — ka. Atunci avem 


A 
(2.7) n(p”) = (p-Da+ Yuri. 

În general pentru 

(2.8) D= Di pap, 

valoarea funcţie 7 este dată de formula: 

(2.9) n(n) = max (n(p*), n(p3*),...,n(ps*)) » 


formulă ce o datorăm lui Kempner [2]. 


Remarca 2.2. Dacă n e N* are descompunerea în produs de numere 
prime (2.8), unde p; sunt numere prime astfel încât p. < p2 < ... < 
ps, iar s > 1, atunci algoritmul lui Kempner de calcul al funcţiei n 
este 


(2.10) p(n) = max ( . (cu) » D2* (aa) (i) i 


„2 Ps? (a ) , 
[ps] (ps) 


unde prin (ap) (2) înţelegem că numărul a este “scris” în baza de nu- 


meraţie generalizată p (notată a|,]) şi este “citit” în baza de numerație 
p (notat 64), unde £ = app)), 6! p. 391. 


3. PROGRAME 


Programul 3.1. Funcţia care numără cifrele în baza de numerație p a 
lui n 
ncb(n,p) := |return ceil(log(n,p)) if n > 1 
return 1 otherwise 

unde funcţia utilitară Mathcad ceil(-) este funcţia parte întreagă su- 
perioară. 
Programul 3.2. Program pentru generarea a bazei de numerație ge- 
neralizată a lui p (notată de Smarandache [p]) pentru un număr cu 
m cifre 
a(p,m) :=|fori e 1..m 

pl 

p=] 

return a 


Qj 4— 
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Programul 3.3. Program de generarea a bazei de numerație p (notată 
de Smarandache (p)) pentru a putea scrie numărul a 
b(a,p) :=|return (U)ifp=1 
for i e 1..ncb(a,p) 
b; — pl 
return b 


Programul 3.4. Programul determină cifrele bazei de numerație ge- 
neralizată |p] pentru numărul n 

Nbg(n,p) :=|m — ncb(n,p) 

a 4 a(p,m) 

retur (1) 4] 10 

fori em..l 


n 
e; — trunc | — 
Qi 


n <— mod (n,a;) 
return c 


Programul 3.5. Program pentru funcţia Smarandache 
n(n) := return ” Err. n nu este întreg” if n Z trunc(n) 
peturi” Er mu La) pl 
return (|)ifn=1 
f<- Fa(n) 
pe fb 
uzi P4z) 
for k = 1..rows(p) 
m <— Pa * Nbg(Qk, Pr) * bla; Dx) 
return max(7) 


Acest program apelează programul F'a(n) de factorizare cu numere 
prime. Programul foloseşte remarca [2.2ja lui Smarandache referitor 
la algoritmul Kempner. 


Programul 3.6. Programul de generarea numerelor prime până la L 
cu ajutorul ciurului lui Atkin. 

Atkin(L) := |eprimr — O 

A — floor(VL) 

for j e L..A 


70 EI 
for k e 1..ceil (272) 


n 4k2 + 3? 
m <— mod(n, 12) 
eprim — =eprima if n < LA(m=1yvm=5) 
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7 RI e 
for k e L..ceil LI) 
n — 3k2 + 3? 


eprima, — —=eprim if n < L A mod(n, 12) = "7 


n) 


for k e j + L..ceil 3 


n 4— 3k2 — 2 
eprim — eprima if n < L A mod(n, 12) = 11 
for k e 5,7,..A 
ori =15, E if eprima 
eprimj.2 4 0 
eprim — 2 
eprima —3 
]—3 
for n € 3,5,..L 
if eprimn 
primej —n 
d del 
return prime 


Programul 3.7. Programul de factorizare a unui număr natural. 

Acest program foloseşte vectorul numerelor prime p care este ge- 
nerat cu ciurul lui Atkin, [1], cel mai rapid program de generare a 
numerelor prime până la limita precizată. Apelul ciurului lui Atkin 
se face cu secvenţa: 


to = time(0) 
p:= Atkin(105) last(p) = 5761455. plasttpy = 99999989 


ta = time(1l) (tb, — to)sec = 82.88sec. 


Fa(m) :=|return ("m =” m” > ca ultimul p*) if m > (prast(p))? 
e aa | 
k +0 
fe) 
while m > pj 
îf mod (m,p;) =0 
k—k+I 


otherwise 


6 CIRA OCTAVIAN 


f — stack|f, (p;, k)] if k > 0 

A a ai 

k —0 
f — stack|f, (p;, k)] if k > 0 
return submatriz(f,2, rows(f), 1,2) 


Programul 3.8. Program de calcul al valorilor funcţiei pentru nume- 
rele naturale până la 10. 
ValFS(N):=|VFS, 1 
for n e 2..N 
VrFS, — n(n) 
return VFS 


Apelul acestui program va genera un vector care va fi salvat ca un 
fişier de tip prn. Astfel vom putea folosi acest fişier pentru rezolvarea 
ecuaţiilor diofantice. Cu ajutorul funcţie Mathcad time s-a contorizat 
timpul de generare al vectorului V/F'S. 


N := 105 to = time(0) VFS:= ValFS(N) ti : time(1) 
(ti — to): sec = "1: 7: 32.623” - hhmmss 
WRITEPRNCVFSprn”):= VFS last(VFS) = 1000000 
Pentru informare se afişează ultimul indice al vectorului V /'S. 


4. REZOLVAREA ECUAȚIILOR DIOFANTICE 


4.1. Ecuația (2069). Pentru căutarea soluţiilor vom citi fişierul V F'S.prn 
cu ajutorul funcţiei Mathcad READPRN 


VFS:= READPRN(..NVFSprn”) last(V F'S) = 1000000 


unde comanda last(V SF) ne indică ultimul indice al vectorului VF'S. 

Generarea fişierului V/F'S odată şi citirea fişierului generat în do- 
cumente Mathcad care determină soluţii ale ecuaţiilor diofantice face 
ca să rezulte o mare economie de timp de execuţie a programului de 
căutare a soluţiilor. 


Programul 4.1. Programul de determinare a soluţiilor ecuaţiei (2069). 
TE 0 2060 (apesi dea Ora Oe da) 9 (ue “ru as) 
u <— last(VFS) 
TOT: duba 
for n € an..ba 
for z € az..ba 
N mM-LĂ+N 
| SA VES = 
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|S — stack|S, (mn a)] if q 
return S 


Apelul programului se face prin secvenţa: 
Bai a ss Ai oi 0 d, > bsi 10 
Numărul total de cazuri verificate este: 
(Drm — Qm + 1)(bn — an + 1)(br — aa + 1) = 1440. 
Apelul programului Ed2069: 
to : time(0) Sol := Ed2069(a, bas Aa ns a,b) bt := time(l) 


Timul de execuţie în secunde şi numărul de soluţii rezultă din: 
(îi —to):s=0011-s  rows(Sol) — 1 = 36 


Pentru m € (2,3,...,10),n e 41,2,...,10Y şi x e 11,2,...,16Y cele 
36 de soluţii ale ecuaţiei diofantice p(m : 2 + n) = z sunt: 


m|n| xi |min| x||m|n| x 
2| 4| 4 4| 7| 7 Za 28 
2| 4| 6 4| 86| 4 7| 5|10 
28-28 4110| 5 7| 6| 6 
2| 5|10 4]10|10 77| 4 
2| 6| 6 5: 44 7[+.8|8 
2| 7| 7 9, 35 8|.:5|.19 
De || ah [ra 5] 6| 6 8| 7| 7 
2|10| 5 5 :Z 838| 9| 9 
31| 5| 5 9-99 9 2 
N [IZA IZA 6| 9| 6 9 |10|10 
3| 7|14 6|10| 5| [10| 7|14 
3| 8| 8 7] 3| 6] [10|10| 5 


Valoarea maximă a soluţiilor x este 15. 


4.2. Ecuația (2070). Pentru căutarea soluţiilor vom citi fişierul V F'S.prn 
cu ajutorul funcţiei Mathcad READPRN 


VFS:= READPRN(..NVFS.prn”) last(V F'S) = 1000000 . 


Programul 4.2. Programul de determinare a soluţiilor ecuaţiei (2070). 
E02070(a+ bras das Oras das ba) > | 5 (m "n a) 
u — last(VFS) 
for m € am..bm 
for n € an..ba 
for z € az..ba 
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N mM-XLI+n 

dp <UAVES = m+n 

S — stack|S, (m n z)] îf q 
return $ 


Apelul programului se face prin secvenţa: 
dpi 2.0 920 "ap bt 20 ag 1 bai 16 
Numărul total de cazuri verificate este: 
(Br — am + I)(bn — am + L)(ba — aa + 1) = 7220. 
Apelul programului Ed2070: 
to : time(0) Sol := Ed2070(am, bms Gas bn; br) bi := time(l) 
Timul de execuţie în secunde şi numărul de soluţii rezultă din: 
(îi —to):s=0.853-s  rows(Sol)— 1 = 14 


Pentru m € (2,3,...,20),n e 41,2,...,20! şi x e 12,3,...,20Y cele 
14 de soluţii ale ecuaţiei diofantice n(m : 2 + n) = z sunt: 
m|2|4|416|6|8|10|12|12|14116|18|18|20 


n|1|1|1|1|1|1| 1| 1| 1| 1| 1| 1| 11|] 
x|412]|614|8|6|12|]10114|12|18|16|20|18 


Valoarea maximă a soluţiilor x este 20. 


4.3. Ecuația (2071). Pentru căutarea soluţiilor vom citi fişierul V F'S.prn 
cu ajutorul funcţiei Mathcad READPRN 


VFS:= READPRN(..NVFS.prn”) last(V F'S) = 1000000 . 


Programul 4.3. Programul de determinare a soluţiilor ecuaţiei (2071). 
EPOTI(asa Des dn 0 daca Bi ) 25 EC m 180.) 
u <— last(VFS) 
for m € am..bn 
for n € an..ba 
for z € az..ba 

N m-rI+n 

dm su AVES,=al 

S «— stack|S, (mn 2) if q 
return $ 


Apelul programului se face prin secvenţa: 
dig te Dia 319 ist ba = 15 Azi SA ba 19 
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Numărul total de cazuri verificate este: 
(Br — am + (ba — aa + (ba — az + 1) = 3990. 
Apelul programului Ed2071: 
to : time(0) Sol := Ed2071(am, bms ns bn a,b) ti := time(l) 
Timul de execuţie în secunde şi numărul de soluţii rezultă din: 
(îi —to):s=0.02-s  rows(Sol) — 1 = 24 


Pentru m € 42,3,...,15),n € 41,2,...,l5işix e 41,2,...,19Y cele 
24 de soluţii ale ecuaţiei diofantice p(m - 2 + n) = z sunt: 


m 2 2 2 3 3 4455 7 8 9 10 
n |= | 3 10 12 7 9 46 1 3 15 129 6 
5 3 3 3333333 3 33 3 


10 11 11 11 12 12 13 13 14 14 15 
15 3 12 15 9 12 6 9 3 6 3 
3 3 3 3 3 3 3 3 3 3 3 


Valoarea maximă a soluţiilor z este 3. 


4.4. Ecuația (2072). Pentru căutarea soluţiilor vom citi fişierul V F'S.prn 
cu ajutorul funcţiei Mathcad READPRN 


VFS := READPRN(..NVFSprn”) last(V F'S) = 1000000 . 


Programul 4.4. Programul de determinare a soluţiilor ecuaţiei (2072). 
Ed2072( ax, br dr br as ba) > |9 ("m Pan 2) 
u — last(VFS) 
for m € am..bn 
for n € an..ba 
for z € az..ba 

4 Ea al big 

d <UAVP Spa" 

S «— stack|S, (mn z)] if q 
return S 


Apelul programului se face prin secvenţa: 
ARI Ag e Pda 92 Dai 0 Mia 32 pt) 
Numărul total de cazuri verificate este: 
(Bm — am + L)(bn — an + L)(bz — ax + 1) = 576. 
Apelul programului Ed2072: 
to : time(0) Sol := Ed2072(am, bas Cn ns Qa; ba) tb := time(1) 
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Timul de execuţie în secunde şi numărul de soluţii rezultă din: 
(îi —to):s=0.1l-s  rows(Sol) — 1 = 12 


Pentru m € (2,3,...,9),n € 42,3,...,9)şiz e 11,2,...,1l0Ycele 12 
de soluţii ale ecuaţiei diofantice p(z”) = x” sunt: 


m 2 3 3455666777 
iese | 3 2 De, 3 DID o Di 2 49 
i 2 2 3 3 5 2 45 2 4 72 


Valoarea maximă a soluţiilor z este 7. 


4.5. Ecuația (2073). Pentru căutarea soluţiilor vom citi fişierul V F'S.prn 
cu ajutorul funcţiei Mathcad READPRN 


VFS:= READPRN(..NVFS.prn”) last(V F'S) = 1000000 . 


Programul 4.5. Programul de determinare a soluţiilor ecuaţiei (2073). 


020710 ba se ba da 0) >) 9 (Cm "0 3) 
u <— last(VFS) 
for m € am..bm 
for n € an..bn 
for z € ag..ba 
q— x" <uA(VFS)” = VPSan 
S «— stack|S, (m n z)] îf q 
return $ 


Apelul programului se face prin secvenţa: 
Bia = 2 Ut 9! dp = bai ap = 2 bai 25 
Numărul total de cazuri verificate este: 
(Bra — Qm + 1)(ba — ap + 1)(be — aa + 1) = 1536. 
Apelul programului Ed2073: 
to : time(0) Sol := Ed2073(am, bras Aa bn Qa; ba) tb := time(l) 
Timul de execuţie în secunde şi numărul de soluţii rezultă din: 


(îi —to):s=0.014-s  rows(Sol) — 1 = 20 
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Pentru m € (2,3,...,9),n e 42,3,...,9) şiz e 12,3,...,25) cele 20 
de soluţii ale ecuaţiei diofantice n(2)” = p(2”) sunt: 


m De 0 De i d Di 2 AD o SD 0 d A 
n |= | 23 3 3 4 4 4 455 5 566 
SE 2 2 36 3 6 8 24 5 8 10 15 45 


Valoarea maximă a soluţiilor x este 24. 
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